Resonating valence bond wave function for the two dimensional fractional spin liquid 
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The unconventional low-lying spin excitations, recently observed in neutron scattering experiments 
on CS2CUCI4, are explained with a spin liquid wave function. The dispersion relation as well as the 
wave vector of the incommensurate spin correlations are well reproduced within a projected BCS 
wave function with gapless and fractionalized spin-1/2 excitations around the nodes of the BCS gap 
function. The proposed wave function is shown to be very accurate for one-dimensional spin-1/2 
systems, and remains similarly accurate in the two-dimensional model corresponding to CS2CUCI4, 
thus representing a good ansatz for describing spin fractionalization in two dimensions. 
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The fractionalization of spin excitations in one- 
dimensional (ID) spin-1/2 antiferromagnetic (AF) 
Heisenberg systems was conjectured by Faddeev and 
Takhtajan more than 20 years ago, and by now the 
concept of "spinous", as spin-1/2 elementary excitations 
for ID quantum AF systems, has been well established 
both theoretically 2] and experimentally JH. 

An important and intriguing issue raised by Ander- 
son after the discovery of high-Tc superconductors 
is whether fractionalized gapless spinous can be defined 
even in higher dimensions. Many theoretical studies have 
been done so far analytically 5] and numerically [^, 
mainly by considering weakly coupled chains. Most stud- 
ies suggest that unconventional ID features are very un- 
likely to occur in higher dimensions since more conven- 
tional states, e.g., magnetically ordered, are stabilized 
with a weak inter-chain coupling. 

The same question can be also addressed experimen- 
tally. A series of recent inelastic neutron scattering ex- 
periments on CS2CUCI4 by Coldea ei 0, showed 
that, as in ID systems |3, the spectrum on this ma- 
terial consists of a broad incoherent continuum at each 
momentum, interpreted as spin fractionalization, despite 
the fact that the system is clearly two-dimensional (2D). 
It was also found that the system is described by the fol- 
lowing 2D spin-1/2 AF Heisenberg model (2DAFHM) on 
the triangular lattice (see Fig.QJa)): 



H = jJ^Si- Si + J' J2 

(iJ> {(iJ» 



(1) 



with the intra-chain coupling J = 0.374 meV and the 
inter-chain coupling J'/ J = 0.33 0- Here the sym- 
bol (i,j) (((i,j))) indicates nearest neighbor sites along 
the chain (between different chains). Additional terms 
such as the Dzyaloshinskii-Moriya interaction and the 
coupling in the third direction, both of the same order 
~ 0.02 meV 0, are certainly relevant to explain the var- 
ious low-temperature phases 0|- However in this study 
we ignore these much smaller terms as our main purpose 
is to determine the minimal model and its ground state 
wave function (WF) which may lead to spin fractional- 



ization in the 2D system. 

In this Letter we propose a projected BCS WF for the 
ground state of Eq. |lj: |$) = :P|BCS), where |BCS) is 
the ground state of the following mean-field BCS Hamil- 
tonian. 
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Here = ek — /i, Ck — — 2cos(k • fi), /i is the chem- 
ical potential, ^ the creation operator of an electron 
with momentum k and spin s = ±1/2, and Ak the 
real gap function for singlet pairing with Ai symmetry. 
The variational state |<I>) for the spin Hamiltonian H is 
obtained by applying to |BCS) the Gutzwiller projec- 
tor V onto the subspace of singly occupied sites. The 
quantities Ak and fi are determined following the varia- 
tional principle, i.e., by minimizing the energy E{\^)) = 
(<i>|i7|<i>)/(<i>|$). For this purpose, the Fourier transform 
Aij of Ak is truncated up to the third distance |i — j| 
along the chain so that 10 independent parameters for 
Ak are optimized along with fj, ^ 10^ . Note that when- 
ever /i is finite, Hqcs no longer possesses particle-hole 
symmetry and therefore |(f>) violates the Marshall sign 
rule 01 . In the following, coupled L-site chains for a 
total number N — L'^ oi sites with periodic boundary 
conditions are used unless otherwise stated. In order 
to explore the quality of |<I>), the variance of the en- 
ergy, ct2(|$)) ^ ($|(i/-£;(|$)))2|$)/($|$), is also cal- 
culated fl^]. For comparison, the same spin liquid WF 
1$) is applied to other models, the uncoupled Heisenberg 
chains with J' — and the 2DAFHM on the square and 
the triangular lattices with J' /J = 0.33 and J'/ J = 1, 
respectively (see Fig.^a)). For each model the gap func- 
tion Ak and fi are optimized independently. 

The main point of this Letter is that the variational 
approach represents an accurate theoretical tool for ex- 
tending the notion of spinous in higher dimensions. With 
this approach we are able to reproduce the qualitative, as 
well as the quantitative features of the experiments, con- 
firming the possibility of spin fractionalization in 2D with 
a transparent theoretical framework. Following Wen , 
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we also show that the excitations of -ffecs are related to 
the physical excitations of the spin Hamiltonian. 
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a similar variational WF with three Bogoliubov modes. 
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FIG. 1: (a) Lattice geometries studied, ri and f2 are the 
primitive translational vectors, (b) Lowest energy spin-1/2 
excitations of the IDAFHM as a function of momentum k. 
Dashed line is the BCS spectrum iSk scaled by a factor to 
match the exact bandwidth, (c) Loci of k-points where = 
(dashed lines) and Ak = (solid lines). The boundary of 
the first Brillouin zone (BZ) for the triangular lattice is also 
denoted by long dashed lines. 



As a first step, let us verify that the highly non triv- 
ial low-energy spectrum of the ID spin-1/2 AF Heisen- 
berg model (IDAFHM) with nearest neighbor coupling 
J is well reproduced by the WF |<i>). It has been well 
known that |$) describes almost exactly the ground state 
WF El [13. Here we show that even the low-lying ex- 
cited states can be constructed from |<i>). For the ID 
model, the optimized fi is zero and therefore the BCS 
spectrum = + Aj^ of Eq. (0) has gapless excita- 
tions at k = ±7r/2, exactly as the spinon spectrum of the 
IDAFHM does 0. Since the elementary excitations of 



Hbcs with energy iJk are described by the Bogoliubov 
modes, 7^ ^, the simplest variational state for the spinon 



at momentum k is |k) = 7^7^. JBCS). To see whether 
this state |k) corresponds to a spinon state, we consider 
a ring with odd number of sites L = 31 and z-component 
of the total spin 5**°* — —1/2. For this case it is known 
that a well defined spinon exists only for half of the total 
Brillouin zone (7r/2 < |k| < tt). As shown in Fig. QJb), for 
this branch, the WF |k) represents fairly well the excited 
state with a spinon at momentum k, as can be verified by 
the good accuracy in energy. Notice that, although the 
projection V is crucial to gain a quantitative agreement 
for the spectrum, the BCS spectrum E]^ already gives a 
qualitatively correct feature of gapless excitations with fi- 
nite spinon velocity at the right momentum k = 7r/2 (see 
Fig^) . It is also possible to obtain an accurate descrip- 
tion for the remaining branch of the spectrum by using 
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FIG. 2: (a) Variance of energy cr^(|$)) iov N = L x L clus- 
ters. UC stands for uncoupled chains, (b) Estimate of the 
variational error in the energy per site. The FN variational 
error for the anisotropic triangular lattice and the variance 
extrapolated error for uncoupled chains are also plotted. 



Even in 2D, the elementary excitations of i?BCS define 
quite naturally, after Gutzwiller projection , th e spin-1/2 
fractionalized excitations of the 2DAFHM [ij. In order 
to consider the possibility of spin fractionalization in 2D, 
it is therefore important to show that the WF |$) can 
appropriately describe the ground state of Eq. ^ . The 
most important effect of frustration on the proposed WF 
1$) is generally to change /z to a nonzero value. This 
should be contrasted to the case of the 2DAFHM on the 
square lattice with no frustration such as the one shown 
in Fig.^a), where the optimized fi is always zero. There- 
fore the Marshall sign rule for the state | $) is violated on 
the triangular lattice, and the gapless excitations of i?BCS 
are no longer located at commensurate points. For exam- 
ple, the optimized value for TV = 18 x 18 is = 0.110(4). 
Thus one can make a spin-1 gapless excitation at wave 
vector Qx = 2cos~^(— /x/2) ~ 1.035(2)7r, which is close to 
the experimentally observed incommensurate wave vec- 
tor on CS2CUCI4 |3- Indeed, as shown in Fig. ^c), the 
nodal lines of Ak intersect the lines defined by = 0, 
and thus four gapless modes exist. 

To check the quality of the proposed WF | $) , the vari- 
ance fT^(|$)) is calculated in Fig. I^l^a), where for com- 
parison the variances are also presented for uncoupled 
chains with J' = (spin liquid) and for the 2DAFHM 
on the square lattice with J'/J — 0.33 (non spin liq- 
uid) . As seen in Fig. EJa) , the variance for the triangular 
case is much smaller than the variance corresponding to 
the square lattice and instead this value is much closer 
to the one for the uncoupled chains. To further explore 
the quality of the WF, we also calculate in Fig. [^b) the 
error of the variational energy from the exact value cq: 
5E ~ i?(|$))/iV — cq. Since the exact energy for the 
triangular lattice case is not known, we estimate eg by 
extrapolating to zero variance the variational energies 
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corresponding to |$) and to the one Lanczos step WF 
l^iLs) = (1 + aH)\^) with optimized a 0]. Indeed this 
procedure works well for the uncoupled chains as shown 
in Fig.|2Ib). Though the WF |<I>) for the triangular lattice 
is not as accurate as for the uncoupled chains, it remains 
of high quality, much better than the square lattice case. 
It should be also emphasized that the corresponding val- 
ues for the best available variational WF with explicit 
AF order of the widely studied isotropic 2DAFIIM on 
the square lattice with J' = J 'Ti| are ct^ ~ 0.015 and 
SE ~ 0.004 in the same unit as in Fig.|21 Also comparing 
these values to ours, the quality of the WF |$) for H is 
exceptionally good for 2DAFIIM's. 

Though the variational approach can provide accu- 
rate upper bounds for the energy, it may suffer the well 
known difficulty that the low energy modes are usu- 
ally insensitive to the total energy, and sometimes com- 
pletely different WF's give more or less the same en- 
ergy [l7|. In order to reduce this difficulty, quantum 
Monte Carlo methods, such as a recently developed fixed 
node (FN) technique, may be used [SlUJi- Within this 
more powerful FN approach, the variational WF |<i>) is 
used to approximate only the signs of the ground state 
WF. In the basis of electron configuration |a;) all the off- 
diagonal matrix elements of the Hamiltonian such that 
Sx'x = ^{x'){x'\H\x)^{x) < are considered exactly. In- 
stead all the ones with Sx'x > are treated semiclassically 
and traced to the diagonal term, defining an effective FN 
Hamiltonian i/pN- Since the semiclassical approxima- 
tion is usually good at low energy, we believe that this 
approach provides a sensible test to examine the stabil- 
ity of the variational WF at low energy. Indeed, if, for 
a properly chosen WF, Sx'x is always non positive for 
X x' , the FN approximation is exact and may change 
completely the low energy properties of the variational 
WF. Namely, for the 2DAFIIM on the square lattice, it 
is possible to obtain the exact AF ordered state provided 
that the variational WF satisfies the Marshell sign rule, 
a condition that does not necessarily imply AF order. In 
the frustrated case no WF's which fulfill Sx'x < for 
X ^ x' are known. Therefore, in order to verify the qual- 
ity of the FN approximation, we estimate by quantum 
Monte Carlo method the fraction w between the posi- 
tive and the negative off-diagonal elements. It is found 
that, due to the quality of our WF, the ratio w remains 
very small (~ 0.03), implying that i/pN contains a large 
fraction (~ 1 — w) of matrix elements equal to the exact 
ones, showing that, in the present case, HpN represents a 
reliable approximation of H. We finally remark that the 
ground state j'l'o) of ifpN can be used as a variational 
state of H with lower energy [li|, £^(|^'o)) < -E(|$)), as 
clearly seen in Fig.El^b). 

In Fig. 121 the spin correlation functions (z-component) 
(^/;|S'^(i)S'^(j)|V')/(V'IV') are calculated for the 2DAFHM 
on both the triangular and square lattices, using the vari- 
ational = 1$)) and the FN = |*o)) mcth- 
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FIG. 3: Spin correlation functions C,;(Z) — {Sz{f)Szir + Iri)) 
as a function of distance I along the fi directions for the 
2DAFHM on the triangular lattice with J' /J = 0.33 (a,b) and 
J' /J — 1.0 (e), and on the square lattice with J' /J = 0.33 
(c,d). Open and solid marks are for the variational (VMC) 
and FN calculations, respectively. Lines are guides to the eye. 



ods. For the square geometry, it is clear that the FN 
approach considerably enhances the large distance corre- 
lations, suggesting that the spin liquid state is unstable 
toward AF long-range order, in agreement with previ- 
ous studies 0] . A completely different scenario is evident 
for the triangular geometry where these long distance 
correlations do not appreciably change with respect to 
the variational |<I>), strongly indicating that the spin liq- 
uid state 1$) is stable, namely close to the exact ground 
state, at least for the clusters studied. 

In order to show further the reliability of our calcu- 
lations, the spin correlation functions for the isotropic 
2DAFIIM on the triangular lattice with J' = J are also 
calculated in Fig. (S^e). Here a "PlBCS) with good varia- 
tional energy is obtained with = 0, and a simple gap 
function with Aij = ±1 restricted up to second near- 
est distances and with properly chosen signs To] . Notice 
that, even by using such a simple spin liquid WF with 
a rather large value of w ~ 22%, the FN method cor- 
rectly reproduces the classical three sublattice antiferro- 
magnetic correlations |20j . 

Encouraged by the above results, we now consider the 
low-lying excitations on the spin liquid state |$) for the 
2DAFHM on the triangular lattice with J' /J = 0.33. To 
this end the FN method is applied. Since w is rather 
small for this case, we expect that this approach should 
give a reliable description of the exact excitation spec- 
trum of H. By using the forward walking technique |2l|, 
one can evaluate the imaginary time (r) evolution of the 
following quantity: 



(^|g,(k)e-"g^^g,(-k)|^o) 
($|e-^-f^FN|^p) 



(3) 
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FIG. 4: (a) A typical semi-log plot of the imaginary time 
correlation function, (b) Calculated spin-1 excitation ener- 
gies (solid squares) along k//ri for the anisotropic triangular 
lattice with iV = 10 x 20, J'/J = 0.33, and J = 0.374 meV. 
(c) The same as in (b) but along k _L ri for different kx- 
Here G2 = A-K/^f%. For comparison, experimental data on 
CS2CUCI4 (open circles) are also shown in (b) 



where 5^(k) = e*'''"S'z(r)/'\/7V. By simple inspec- 
tion, the lowest spin-1 excitation energy E^^^ of i/pN can 
be obtained by fitting the large r behavior of S'(k, t) cx 
e"^"-; . An example of how the method works is given in 
Fig. EI a). For large r, ln5(k, r) is almost linear and thus 
we can estimate the lowest excitation energy for each k. 
The results are summarized in Fig. QJb) and (c). While 
the spectrum has a strong ID characteristic, it shows a 
visible and non trivial dispersion in the k-direction per- 
pendicular to Ti. Moreover, as seen in Fig. EJa), the 
linear fit is not perfect in the intermediate times r J ~ 1 
implying that there are higher energy (~ J) contribu- 
tions to the spectrum, suggesting that the full spectrum 
of the dynamical structure factor may be consistent with 
the incoherent continuum observed in the experiments. 
The extremely good agreement between the calculated 
spectrum and the available experimental data 7'] , shown 
in Fig. ^b) , indicates that the 2D model of Eq. ^ de- 
scribes correctly the experiments, thus supporting the 
appearance of a spin liquid state. Since no sizable mag- 
netic order was found in this 2D model (Fig.OJ, it is very 
likely that the experimentally observed ordered phases 
at very low temperatures should be due only to 3D ef- 
fects j8j. We cannot exclude the possibility that a small 
spin gap could remain in the thermodynamic limit, as 
our numerical resolution is limited in the available finite 
size clusters. 

In conclusion, we have shown that, by studying sev- 
eral models, a simple variational approach is capable of 
describing the ground state properties as well as the low- 
lying excitations of a spin liquid not only for ID systems 
but also for a 2D frustrated spin- 1/2 model. Within this 



approach it appears that, when the i?BCS spectrum re- 
mains gapless at commensurate k as in a non frustrated 
square lattice, the spin fractionalized state undergoes a 
magnetic instability. By contrast, in the frustrated trian- 
gular case the i?BCS spectrum becomes gapless at incom- 
mensurate k as soon as J'/J > 0, and it appears much 
more difficult to destabilize the fractionalized state by a 
magnetic phase transition. It should be emphasized that 
both the experimental and the present numerical work 
consistently suggest that a spin liquid with gapless spin- 
1/2 excitations appears possible in 2D, its stability being 
intimately related to the incommensurate momenta of 
the gapless modes. 

Our results may have some impact even for high-Tc su- 
perconductors, where the carrier doping naturally leads 
to incommensurate nodal Fermi points. The existence 
of a nodal spin liquid has been also conjectured before 
for interacting Bose systems |22b i23j where a Bose metal 
is the natural bosonic counterpart of a spin liquid with 
gapless spin excitations and no magnetic order. 
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by INFM-MALODI and MIUR-COFIN 2001. 
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